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Abstract 

In a previous work one of the authors proposed a simple model for study- 
ing systems under pressure based on the Thomas-Fermi (TF) model of single 
atom. In this work we intend to extend the previous work to more gen- 
eral Thomas-Fermi models where electronic exchange and correlation are in- 
troduced. To do so, we first study numerically the equation obtained by 
H.W.Lewis (TFDL) which introduces the effects of exchange and correlation 
into the original TF equation; next the procedure followed in the previous 
work is extended to the new approach and a specific example is illustrated. 
Although one could expect that no big differences were produced by the gen- 
eralized TF model, we show the qualitative as well as quantitative equivalence 
with detailed numerical results. These results support the robustness of our 
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conclusions with regards to the model proposed in the previous work and 
give the character of universality (i.e. to pass from one atom to another, 
the quantities calculated must be simply scaled by a numerical factor) to the 
properties of compressed systems shown in this work. 
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I. INTRODUCTION 



The statistical model of atom developed by E.Fermi [1] and known as Thomas-Fermi 
model, although based on a highly simplified theoretical framework, has been proven sur- 
prisingly good in predicting basic properties of condensed matter [2]. A particular feature 
of such a model is the description of compressed atoms, in this case the predicted properties 
were confirmed experimentally [3]. In a preceding publication [4], we focused the attention 
on this particular aspect of the theory and proposed a simple model for describing systems 
under pressure. The central point of that work was the development of the concept of 
"statistical ionization". In simple terms, this is a universal semianalytical function which 
enables one to describe, as a function of the distance from the point-like nucleus, the balanc- 
ing process between the antibinding and binding contribution to the total energy within the 
compressed atom. In spite of the approximations done and extensively discussed, we under- 
lined the utility of the proposed model as a tool for investigating at a basic level and at low 
computational cost, properties of systems under pressure. However it was also underlined 
that the properties of semianalycity and universality of the function disappear when higher 
order of approximation are introduced in the basic TF theory. In the light of what stated 
before, in this publication we intend to extend the treatment of the previous work to more 
sophisticated models of the TF approach. That is to say, to include the effects of exchange 
and correlation into the original TF model and obtain in this framework a "generalized sta- 
tistical ionization function". The extension of the original TF model is due to P.A.M.Dirac 
(TFD) [5] for the exchange while for the part relative to the correlation, several approaches 
has been proposed; in this work the one proposed by H.Lewis (TFDL) [6] has been chosen 
because it is simple and appropriate for the compressed case, since the treatment of the elec- 
trons becomes exact in high density limit (see appendix). The paper is organized as follows; 
first we obtain the TF, TFD and TFDL equation within a single generalized approach, then 
we show numerical solutions for different atomic numbers in the neutral uncompressed case, 
next we illustrate numerical results for the "generalized statistical ionization function" in 
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case of Z = 50 {Z atomic number). Finally the equation of state of a compressed system is 
calculated. Comments on the results obtained as well as on advantages and limitations of 
this model,conclude the work. 

II. A STATISTICAL APPROACH: THE THOMAS-FERMI-DIRAC-LEWIS 

MODEL 

In this section we derive a unique general form for the TF, TFD and TFDL equation. 
In a semiclassical approximation the local electronic density of states can be defined as (see 
for example [7]): 

dn\p{r),{r)] = ^p'{r)dp(r). (2.1) 

The local electron density is therefore given by: 

n(r) = I ^/(r)dp(r) = (2.2) 

which determines the local Fermi vector Pf(i')- 

In the spirit of the (generalized) Thomas-Fermi approach we can express the one-particle 
energy E\p{r),r] through an effective Hamiltonian i7ei[p(r)] of free (interacting) electrons 
plus an electrostatic field V{r) arising from the direct electron-nucleus and electron-electron 
interaction (we must remind that the nucleus is considered a positive point-like charge) : 

E\p{r),r]^HM^)]-eV{r), (2.3) 

For a system at the equilibrium the chemical potential /^(r), defined by the maximum energy 
of max{E[p{r), (r)]} = E[pF{r),r], has to be a constant independent of r. Putting for 
convenience jj, = —cVq we obtain therefore: 

H,,{pF{r)]=e[V{r)-Vo]. (2.4) 

The effective field V^(r) will be determined self-consistently via electrostatic considera- 
tions by using the Poisson's equation: 
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1 n(r) 
V'[V{r) - V,] = {r [V{r) - V^]} = e^, (2.5) 

where r = |r|. Eq. (2.5) can be quite simplified by introducing the variables 0(r) and v{r) 
defined by: 

Ze 



-Hr) = V{r) - Vo, (2.6) 



v'^{r) = {2Ze'^m)-\4:7ieo)rpl{r). (2.7) 

Eq. (2.5) becomes thus: 

0"(r) = ^ (^) i2ZE'mf'iA^eor/'v%r)r-'/^ (2.8) 



and finally, by scaling r — > Paox, we have: 

0"(a;) = v^{x)x-^/^, (2.9) 
this is the "generalized TF equation" , where 



«o = tVS' (2-10) 



is the atomic Bohr radius and 



A closed system of differential equations is now obtained by expressing back the variable 
(j){x) as function of v{x). To this aim the explicit expression of the electronic Hamiltonian 
Hei[p{r)] is needed. Of course this task is in principle quite hard since the solution of a many- 
body system is necessary, and some kind of approximation is required. The original TF 
model represents the simplest approximation, i.e. the Hamiltonian H\p{r)] is approximated 
with the only kinetic term: 

Hl'lpir)] - (2.12) 
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The total Hamiltonian contains therefore only kinetic and the direct Coulomb electron- 
electron interaction, totally neglecting any quantum contribution. By using again the re- 
duced variables 4>{x) and v{x), we can write Eq. (2.12) in the compact form: 

(p{x) = v^{x), (2.13) 

which, together with Eq. (2.9), defines the statistical solution of the Thomas- Fermi model. 

AS stated before higher order of approximation have been introduced in literature to 
include quantum exchange and correlation contribution to Eq. (2.12). Most famous is 
the so called Thomas- Fermi- Dirac (TFD) model which explicitly takes into account the 
exchange energy (see [5,2]). The inclusion of the correlation term has been a more debated 
issue, and different approaches have been proposed. A particular suitable and simple one 
is the interpolation formula given by H.W. Lewis, which becomes exact in the both the 
high-density and low-density limits [6]. We referred and will refer to it as the TFDL model; 
in this case the Hamiltonian contains kinetic, exchange and correlation terms: 



H-TFDLr / _ r(0 e^p(r) me^(l^hi2) 



1 + 



(O.SQttTT — l)np{r)ao 
(1 -ln2)^ 



(2.14) 
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a=[-) . (2.15) 



where 

.9ttJ 

Introducing the variables 0(x) and v{x), Eq. (2.14) reads: 

(f){x) — v'^{x) — av{x)x^^^ — pa^xln aa + v{x)x~^^'^ , (2-16) 

where 



l-ln2 , , 

^^2(0.89a7r-l)' ^^"^^^ 

and 



p = -(l-ln2). (2.19) 

Eq. (2.16) can be considered a generalization of the TF model [Eq. (2.13)] to include 
exchange and correlation energies. We find Lewis's formula particularly appealing since we 
are mainly interested in the high-density regime of compressed atoms where the interpolation 
formula provided by Eq. (2.14) is expected to work quite well. In spite of the simplicity of 
the model, to our knowledge, no explicit numerical or analytical treatment of this equation 
has been published, although several authors [8-11] cite this approach in their work. The 
main reason why this happens is due to the fact that since the TFD model does not produce 
evident changes with respect to the simple TF solutions, the TFDL is expected to not be 
much different. 

However in our paper we derive for the first time explicit numerical solutions of the TFDL 
model and we use it as basic starting point of our analysis; the results will be compared with 
the simple TF model. Although we do not expect large differences with respect to the TF, 
we carry on this study to confirm a qualitative and quantitative equivalence between TF 
and TFDL and later on the consequent advantages of such a conclusion will be underlined. 

III. NUMERICAL STUDY OF TFDL EQUATION AND "GENERALIZED 

IONIZATION FUNCTION" 

In this section we illustrate some TF,TFD and TFDL solutions for the neutral uncom- 
pressed case, the trend shown for this case is preserved for the compressed case. Then we 
show within a generalized TF approach, the procedure to obtain the "statistical ionization 
function" via the total energy as a function of the distance from the nucleus. Eqs. (2.9) and 
(2.16), or Eqs. (2.9) and (2.13), represent two sets of coupled differential equations which 
can be solved to obtain (f){x) or, equivalently, v{x), and as a consequence the local electron 
density n(r). The boundary conditions required to determine a unique solution are provided 
by the asymptotic limit: 
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limc!)(x) = 1 



(3.1) 



and by fixing the radius of a neutral atom: 



0'(xo) = Xo0(xo). 



(3.2) 



Altliougli it lias been widely pointed out in literature, it is important to remind that the 
inclusion of exchange and correlation terms breaks down the universality of the solution. 
Namely, it is not possible to obtain local electron density n(r) for a given atomic number 
Z just by a proper scaling of the physical quantities (length, energy etc . . . ). In Fig. 1 we 
show the radial distribution of 0(x) for a neutral uncompressed atom for the Thomas-Fermi 
model and for the Thomas-Fermi-Dirac-Lewis model with statistically relevant values of 
atomic number: Z = 10, 50, 80. Only a slight discrepancy is found as result of the inclusion 
of the exchange and correlation terms, pointing out the almost quantitative accuracy of the 
TF model. As expected, such a discrepancy vanishes as the atomic number Z is increased. 
This result show that TF and TFDL are quantitatively as well as qualitatively equivalent 
(the trend as stated before holds in the compressed case as well) . 

By the knowledge of the solution ^(a;) we can now evaluate the semiclassical expression of 
the local density of energy E{r) which is related to the one electron energy E\p{r), r] through 
the relation E{r) — J dn\p{r), (r)]£'[p(r), r]. We can identify three different contributions: 



E{r) = E^^{r) + E,_i{r) + E,_,{r), 



(3.3) 



where 




(3.4) 



^e-n(r) 



Ze^ 
47reor 



n(r). 



(3.5) 




(3.6) 



where of course 
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n(r) = ldn\p{r),{r)]. (3.7) 

A particularly meaningful quantity is also the integrated energy defined E{R) by: 

E(R)^ [^inr^drE(r), (3.8) 
Jo 

where i? is a generic distance R < tq. E{R) represents therefore the total energy contained 
in the atomic volume enclosed within the distance R from the nucleus. 

Just as E{r), we can think E{R) as composed by a kinetic, a Coulomb electron-nucleus 
and a Coulomb electron-electron term. After some straightforward manipulations (see also 
[12]) , by using once more the reduced variables, we can write each of these contributions 
as: 



OTT 



me 



me 



X 



dxv^{x)x' 



-1/2 



X 



dxv^{x)x ^^"^^ 



Ee-e{X) 



37r 



me 



2(47reo^)^ 



X 



(3.9) 



(3.10) 



dxv'{x)x'/^ \- r v\x')x'^''' + rv\x')x'-''^ 
ix Jo Jo 

(3.11) 



where X is the scaled variable related to R through the relation it! = PaoX and where we 

made use of the relation: 

/(r)/(r') , , , ^« . . . n 



f f /iW^^^rdr'^ / drATrr^f(r)\- r dr'A7rr'^f(r')+ r dr'Airr'fir') 

J\r\<RJ\r\<ro \r - r'\ Jo ^ ^ W Jo ^ ' Jo ^ ' 



(3.12) 



valid for an isotropic /(r). 

As it was done in [4], we consider R a variable distance within the atom, and E{R), 
as stated before, is the total energy in classical terms, of the sphere of radius R inside the 
spherical atoms. In simple terms, wc do divide the atom into infinitesimally thin concentric 
shells; the energy at R is the sum of the contributions of all of the shells inside R. As a 



consequence the distance at which E{R) has its minimum can be interpreted as the distance 
at which the binding (nucleus-electrons) and antibinding (electron-electron and kinetic en- 
ergy) contributions to the total energy are in exact balance. This allows us to define a sort 
of electron ionization criterion, where the term "ionization" stands for electrons with zero 
or positive energy; thus we can address to E(X) as the "generahzed statistical ionization 
function". The exact meaning of "statistical ionization" with the related limitations has 
been extensively discussed in our previous work, for such a reason we do not spend more 
discussion about it and focus the attention onto the numerical results. 
The number of ionized electron is therefore given by 



Differently from the procedure adopted in our previous work, in this case E{X) can be 
studied only numerically, since it is not possible to solve analytically the integrals and so 
obtain a semianalytical expression. Moreover the dependence of (l){x) on Z does not allow 
to obtain a universal function where to pass from one atom to another is possible by simply 
scaling in Z. Next we show results obtained by studying numerically the "generalized 
statistical ionization function" for the case Z — 50; however since we did find quantitative 
agreement between TF and TFDL, the results shown are valid for any atomic number by 
opportunely scaling them. 



In this section we show results obtained for the "generalized statistical ionization func- 
tion" for the case Z — 50. As it was already underlined, we expect that the results are 
quantitatively and quahtatively equivalent in case TF or TFDL model is used, since the 




(3.13) 



which, written in reduced variables, reads: 




(3.14) 



IV. RESULTS 
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solution of the TF and TFDL equation do not differ, as it was shown in the previous sec- 
tion. This means that the "generahzed ionization function" is the numerical equivalent of 
the semianalytical and universal "statistical ionization function" of our previous work; thus 
the results (and the equation of state) shown here for Z — 50, can be easily generalized to 
any Z by simply scaling. Indeed the function E{X) does not show a different behavior when 
TF or TFDL are used; in figure 2 the plots obtained using TF and using TFDL coincide. 
As it is possible to see in figure 2 a minimum is always obtained for any compression, and 
in figure 3 a particular compression has been chosen and the determination on the number 
of ionized electrons is pictorially illustrated. 

Accordingly to the procedure of our previous work, at this point we can therefore model 
the compressed atom as a core with radius i?core plus a number of ionized electrons A^ion 
which in principle is free and is spread over all the atomic volume. The pressure is therefore 
determined by the density of ' 'free" (ionized) electrons as we proposed in the previous work 
(again limitations are extensively discussed there): 



where Aion is given by Eq. (3.14) and the atomic volume is simply Kitom = t\. First we 
study the ionization as a function of the compression for the single atom (see Fig.4) and 
then we extend the procedure to a macroscopic level so that the resulting equation of state 
is shown in Fig. 5, where the pressure in units of Rydberg over the Bohr atomic volume 
l^ohr = is plotted as function of the volume. 



We generalized the concept of "statistical ionization function" for electrons in a com- 
pressed atom, obtained in a previous work, within the TF approach. The generalization 
was developed by extending the previous analysis to more sophisticated TF models where 
exchange and correlation are considered. We found that there is no qualitative and quan- 
titative difference between the original and the generalized approach. Although one could 




(4.1) 



V. DISCUSSION AND CONCLUSION 
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have expected this result, we exphcitly proved such an equivalence, and for the first time, 
solutions of the TFDL equation where shown. The important conclusion of this work con- 
cerns the fact that because of the equivalence shown, the results reported for the particular 
case of Z = 50 are valid for any atomic number, provided that a rescaling ^^^j where Z* is 
the atomic number wished, is applied; above all the results legitimate the use of the semi- 
analytical and universal ionization function obtained in our previous work with the evident 
advantages of the low computational cost and its extreme simplicity and immediacy. More- 
over the fact that Lewis' formula is exact for the high density limit, means that in our case 
the effects of correlation were well described, as a consequence one can conclude that such 
effects are not relevant for describing atoms under pressure, at least in first approximation; 
this result it is not obvious in an a 'priori analysis. As it was discussed in the previous work, 
the equation of state is a simplification and as a consequence far from being rigorous; for 
instance an open problem of our model is the distribution of the "ionized" electrons which in 
the present work is considered simply uniform. However, due to its simplicity and feasibility, 
the model can be applied for a basic study of compressed matter not only via the deter- 
mination of the equation of state, but also as a basis for developing more efficient analysis 
within more sophisticated theories. The example shown in our previous publication, was 
the determination of the "ionized" electrons for a certain compression and the consequent 
description, in an ah initio method, of the ionized electrons via plane-waves wavefunctions 
while the other electrons can be represented as a core or as localized orbitals; this would 
speed up the convergence for self-consistent calculations of compressed matter. Here we 
can say more, as it is well known, the ah initio calculations are based on the concept of 
pseudopotential; only the valence electrons are explicitly taken into account while the rest 
are placed in a core described by an opportune potential which interacts with the valence 
electrons. Valence electrons and core are known for the uncompressed atoms but one may 
ask what happens in case the system is under high pressure, in this case the valence electrons 
and the core should be redefined according to the degree of compression. In this case our 
model, at basic level' would be very helpful for estimating the number of ionized (valence) 

12 



electrons and define what is the core. The evident advantage of such a procedure stays in 
the simphcity of such an estimate and in the very low computational cost. In conclusion we 
think that the study performed in this work furnishes important information and tools for 
a computational inexpensive and well founded basic analysis of compressed matter. 



VI. APPENDIX 

A. Correlation Energy within Lewis' approach 

We briefly illustrate the procedure followed by lewis to develop a suitable formula for the 
electron correlation, for more details see [6] and references therein. Consider an electronic 
Fermi gas whose maximum momentum is the Fermi momentum pp and density n = 
The correlation energy at high density can be calculated via the Gell-Mann's scheme and 
leads to the following expression 

Ecorr — — (me^/7r^^^)(l — ln2)ln{pF) + const (6-1) 

m and e are respectively the electron mass and the electron charge. This expression must 
be modified in such a way that its validity could be reasonably extended to any density and 
in particular must be exact at low density. At this point Lewis invokes Wigner procedure 
for the calculation of electron correlation for a dilute gas; i.e. one simply needs to note that 
a very dilute electron gas in the ground state crystallize into a body-centered cubic lattice 
and at this point the correlation energy can be determined exactly via a Madelung type 
technique. The expression for low density obtained is 

Ecorr = -(O.SQaTT - l)e'^pF/7rh (6.2) 

where a — (4/97r)^/^. Finally the procedure leads to what Lewis defines as "a suitable 
interpolation formula" for the correlation energy valid for any density: 

(0.89Q;7r — l)pi?7rao 



__me^(l-/n2) 



(1 - ln2)h 
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(6.3) 



where Qq — fi? /me^. This is the expression used by Lewis to incorporate the correlation 
effects into the TFD model. 
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FIGURES 

FIG. 1. Radial distribution of as function of x for a neutral uncompressed atom. Dashed 
line represents the Thomas-Fermi model, solid lines the Thomas- Fermi-Dirac-Lewis one with atomic 
number (from lower to upper line) Z = 10,50,80. The Z = 50 and Z = 80 cases are strongly 
overlapping and are hardy distinguishable in the plot. 

FIG. 2. Integrated energy E{R) as function of R for the Thomas-Fermi-Dirac-Lewis model 
{Z = 50) and different degrees of compression parameterized by the atomic radius (from lower to 
upper line): xq = 1.5, 1.4, 1.3, . . . ,0.9,0.8. E{R) always presents a minimum at a certain R. 

FIG. 3. Top panel: particular of Fig. 2 for xq = 1.2. The minimum of E(R) determines the 
"core" radius -Rcore- Bottom panel: graphical representation of the number of ionized electrons 
Nion corresponding to the above case. Nion is determined by integrating the electron density from 
the "core" radius to the atomic boundary. 

FIG. 4. Number of ionized electrons (solid line, left scale) and core radius Xcore (dashed line, 
right scale) as function of the atomic radius xq. 

FIG. 5. Pressure vs. volume dependence for the statistical ionization model. 
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